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We report a limit on the ultra-high-energy neutrino flux based on a non-detection of radio pulses 
from neutrino-initiated particle cascades in the Moon, in observations with the Parkes radio telescope 
undertaken as part of the LUNASKA project. Due to the improved sensitivity of these observations, 
which had an effective duration of 127 hours and a frequency range of 1.2-1.5 GHz, this limit extends 
to lower neutrino energies than those from previous lunar radio experiments, with a detection 
threshold below 10 20 eV. The calculation of our limit allows for the possibility of lunar-origin pulses 
being misidentified as local radio interference, and includes the effect of small-scale lunar surface 
roughness. The targeting strategy of the observations also allows us to place a directional limit on 
the neutrino flux from the nearby radio galaxy Centaurus A. 

PACS numbers: 98.70.Sa 


I. INTRODUCTION 

The spectrum of cosmic rays (CRs) extends into the 
ultra-high-energy (UHE; > 10 18 eV) regime, in which 
their origin is currently unknown. Above a threshold en¬ 
ergy of ~ 6 x 10 19 eV (for CR protons) they lose energy 
through interactions with the cosmic microwave back¬ 
ground [1, 2], known as the Greisen-Zatsepin-Kuzmin 
(GZK) effect. This results in a dramatic steepening in 
the CR spectrum above this energy [3]. The few CRs ob¬ 
served above the GZK threshold energy must originate 
within the local universe, within a few times the energy- 
loss distance (tens of Mpc). 

CRs across a range of energies are known to be pro¬ 
duced in regions of magnetic turbulence such as super¬ 
nova remnants [4], where charged particles can be ac¬ 
celerated through interactions with magnetic fields. Be¬ 
cause this process starts with a low-energy particle and 
increases its energy to produce a CR, models of this type 
are described as bottom-up models, with the details de¬ 
pending on the particular acceleration mechanism and 
the environment of the source object. In the UHE regime, 
it is proposed that CRs could also originate via a top- 
down mechanism in which UHECRs are decay products 
from hypothetical superheavy particles, with the details 
of the models depending on the assumed properties of 
these particles. 

Different UHECR origin models predict a range of dif¬ 
ferent UHECR source spectra; however, the observed 
spectrum is modified by energy-dependent GZK atten¬ 
uation and photo-disintegration of CR nuclei, making 
discrimination difficult. Different models also predict dif¬ 
ferent astronomical objects from which UHECRs should 
originate, but charged UHECRs are deflected in cos¬ 
mic magnetic fields, obscuring their directions of origin. 


These deflections are less significant for more energetic 
UHECRs, and there is a correlation between UHECR ar¬ 
rival directions and the positions of nearby active galac¬ 
tic nuclei (AGN), a commonly-proposed class of source; 
but these are in turn correlated with the distribution of 
matter in the local universe, and thus with sources for 
other models [5]. There is also a tentative correlation of 
UHECR arrival directions on large scales (~ 20°) with 
the nearby AGN Centaurus A [6]. 

Another approach for exploring the origin of UHECRs 
is to search for counterpart neutrinos at similar energies. 
Both bottom-up and top-down origin models predict an 
associated UHE neutrino flux. In addition, independent 
of the origin model, UHECRs undergoing GZK interac¬ 
tions produce pions which decay to produce additional 
UHE neutrinos. Since neutrinos do not interact signifi¬ 
cantly with intervening matter or radiation fields, these 
neutrinos retain information about the source UHECR 
spectrum [7]. In addition, since they are uncharged, they 
are not deflected by magnetic fields, and their arrival di¬ 
rections correspond directly to their directions of origin. 
Consequently, the detection of even a single UHE neu¬ 
trino (with sufficient angular resolution) could identify a 
source of UHECRs, and sufficient detections to measure 
the UHE neutrino spectrum would provide a strong test 
of competing UHECR origin models. 

Efforts to detect UHE neutrinos focus on searching for 
nanosecond-scale radio pulses produced when they in¬ 
teract in dense media: the interaction initiates a particle 
cascade which develops an excess negative charge primar¬ 
ily through entrainment of electrons from the medium, 
resulting in coherent radiation at wavelengths larger than 
the width of the cascade (A > 10 cm; i.e. radio), directed 
at the Cherenkov angle. This burst of coherent radiation 
was predicted by Askaryan [8], and we refer to it as an 
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Askaryan pulse. The effect requires that the particle cas¬ 
cade occur in a radio-transparent medium; it has been 
observed in laboratory experiments with target media of 
silica sand [9], rock salt [10] and ice [11]. 

UHE neutrino radio-detection experiments have been 
conducted with terrestrial ice as a target medium, mon¬ 
itored with radio antennas embedded within the ice 
[12, RICE], suspended from a high-altitude balloon [13, 
ANITA], or mounted on a satellite [14, FORTE]. They 
have also been conducted with the target medium be¬ 
ing the lunar regolith (with properties similar to silica 
sand), monitored with ground-based radio telescopes as 
suggested by Dagkesamanskii and Zheleznykh [15]. Due 
to the size of the Moon, this latter technique offers the 
largest potential experimental aperture; however, it re¬ 
quires the Askaryan pulse to be bright enough to be vis¬ 
ible at the distance of the Moon (~ 3.8 x 10 8 m), which 
makes it sensitive only to the most energetic neutrinos. 
This makes it well-suited to testing top-down UHECR 
origin models, as these typically predict neutrino spectra 
extending to higher energies than those from bottom- 
up models or GZK interactions [e.g. 16-18]. In addition 
to neutrinos, this technique has some capacity to detect 
UHECRs directly [19, 20], but in this article we consider 
the sensitivity to neutrinos only. 

The LUNASKA (Lunar UHE Neutrino Astrophysics 
with the Square Kilometre Array) project aims to de¬ 
velop this technique for future use with the Square Kilo¬ 
metre Array [21], a radio telescope with a planned sensi¬ 
tivity exceeding all present instruments. As part of this 
project, we have conducted an experiment employing this 
technique with the Parkes radio telescope, searching for 
lunar-origin Askaryan pulses as a tracer of UHE parti¬ 
cle cascades in the lunar regolith. Sec. II contains an 
overview of the design of our experiment, and Sec. Ill 
summarizes our observations and analysis; full details of 
our experimental procedure are published in a separate 
article [22]. As in previous experiments of this type, we 
do not detect any lunar-origin radio pulses, which es¬ 
tablishes a limit on the UHE neutrino flux, through the 
procedure described in Sec. IV. Due to our observational 
strategy, we are able to place limits both on the diffuse 
flux (Sec. V) and on the directional flux from Centau- 
rus A (Sec. VI). Finally, Sec. VII summarizes our major 
results. 


II. EXPERIMENT DESIGN 

Our experiment was conducted the Parkes radio tele¬ 
scope, a single 64 m parabolic antenna, using its 21 cm 
multibeam receiver [23], with a radio frequency range 
of 1.2—1.5 GHz. This receiver is capable of operating 
thirteen simultaneous beams on the sky, of which we 
used four beams at any one time, as shown in Fig. 1. 
Two or three of these beams were placed near the limb 
of the Moon, matching the expected position of a de¬ 
tectable Askaryan pulse [24], while the remaining one or 



FIG. 1. Typical pointings (solid shading) of telescope beams 
relative to the Moon in our observations, with approximate 
system temperatures T sys for each beam. The limb beams are 
more sensitive than the half-limb beam; the off-Moon beam 
is even more sensitive, but useful only for identifying RFI. 
Crosses in each beam indicate the orientation of the orthog¬ 
onal linear polarizations. The receiver orientation was main¬ 
tained such that Centaurus A was always in the direction 
shown. For some of the observing time, the half-limb beam 
was replaced with an additional off-Moon beam (striped). 


two beams were directed away from the Moon and used 
to identify spurious signals from radio frequency inter¬ 
ference (RFI). One limb beam was always maintained at 
the point on the limb closest to Centaurus A in order 
to maximize the sensitivity to neutrinos (or cosmic rays) 
from that direction, as done by James et al. [25]. 


Real-time processing of the signal, after downconver- 
sion to an intermediate frequency range of 50-350 MHz, 
was carried out with the Bedlam backend [26]. A buffer 
of 4 /us of sampled baseband data was maintained for 
each of the eight input channels, corresponding to two 
polarizations on each of four beams. If a peak exceeding 
a trigger threshold was detected in any on-Moon beam, 
and no peak was simultaneously detected on any other 
beam, the contents of these buffers were copied to per¬ 
manent storage. This latter anticoincidence condition 
excluded the majority of RFI pulses, which typically en¬ 
tered from sidelobes of the antenna, and so were detected 
on multiple beams. 
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A. Signal optimization 

The amplitude of a lunar-origin Askaryan pulse is de¬ 
creased by dispersion as it passes through the ionosphere. 
The degree of this dispersion is determined by the iono¬ 
spheric total electron content (TEC); specifically, by the 
column density of plasma along the line of sight to the 
Moon, called the slant TEC (STEC) and typically mea¬ 
sured in TEC units (TECU). The Bedlam backend imple¬ 
ments a digital dedispersion filter which compensates for 
this, but was limited by our imperfect knowledge of the 
STEC. The maximum deviation between the filter set¬ 
ting, based on real-time ionosonde measurements, and 
the more precise STEC measurements available retro¬ 
spectively from Global Positioning System (GPS) data, 
was 9.0 TECU, which determines the real-time loss of sig¬ 
nal strength from dispersion. The loss of signal strength 
in full retrospective processing is determined by the max¬ 
imum uncertainty in the GPS-derived STEC, which was 
3.8 TECU. The maximum value of the STEC at any point 
during our observations was 23.5 TECU. 

Digital sampling at a finite rate can reduce the am¬ 
plitude of an analog pulse, as the sampled points will 
not generally coincide with the peak of the pulse [25]. 
This effect can be counteracted by interpolating between 
the sampled points. The Bedlam backend implements 
two-fold interpolation, reconstructing one point between 
every pair of sampled points, effectively doubling the real¬ 
time sampling rate. In retrospective processing, we per¬ 
formed 32-fold interpolation, almost completely negating 
this effect. 

The phase of a pulse is randomized by the process of 
frequency downconversion, with a consequent potential 
loss of amplitude [26]. In real time, no effort was made 
to compensate for this effect. In retrospective processing, 
we formed the signal envelope, which restores the full 
amplitude of the pulse for any possible phase, at the cost 
of increasing the noise level. 

The individual contributions of the effects described 
above to the signal loss are shown in Table I, for both 
real-time processing and retrospective processing, com¬ 
pared to the case if no processing had been applied. The 
real-time processing was sufficient to cause any signifi¬ 
cant event to exceed the trigger threshold, allowing it to 
be recovered for later analysis. The retrospective pro¬ 
cessing applied to these stored events was sufficient to 
reconstruct almost the maximum pulse amplitude, with 
a maximum signal loss of only 0.4%. 

III. OBSERVATIONS AND ANALYSIS 

Our observations were conducted in runs of 3-5 days 
once per lunar orbit, when the Moon was closest to Cen- 
taurus A, with a total of six runs from April to September 
2010. After subtracting time spent calibrating or other¬ 
wise excluded, we observed the Moon for a total of 148.7 
hours. We recorded candidate events typically at a rate 


TABLE I. Simulated loss of signal strength due to ionospheric 
dispersion, finite sampling rate and unknown pulse phase; 
with no processing applied, with the real-time processing used 
in this experiment, and with the retrospective processing ap¬ 
plied to recorded data. Dispersion is based on the listed 
STEC; see text for details. 


STEC Maximum signal loss 


Processing 

(TECU) dispersion sampling 

phase 

total 

none 

23.5 

15.0% 

21.6% 

17.9% 

41.9% 

real-time 

9.0 

2.3% 

5.6% 

17.9% 

23.1% 

retrospective 

3.8 

0.4% 

< 0.1% 

< 0.1% 

0.4% 


of 1-2 Hz, with a total of 794,568 events stored for later 
analysis. The majority of these events contained only 
fluctuations from the thermal noise background, but the 
events containing high-significance peaks were dominated 
by RFI pulses, despite the majority of such pulses being 
excluded by our real-time anticoincidence filter, at a ra¬ 
tio of ~ 150:1. For full details of our observations, and 
of the analysis described below, see Ref. [22] 

After applying the full retrospective processing de¬ 
scribed in Sec. IIA to optimize the signal-to-noise ratio, 
we applied a series of cuts to remove the RFI events. 
These retrospective cuts refined the real-time anticoinci¬ 
dence filter, applying more stringent criteria to exclude 
pulses detected in multiple beams, at the cost of in¬ 
curring an increased probability of incorrectly excluding 
non-RFI events. We also applied cuts to exclude events 
with strong narrow-band RFI, which could saturate the 
analog-to-digital converters; events with a pulse width in 
excess of 10 ns, which is significantly greater than the 
expected duration of a band-limited Askaryan pulse; and 
events which occurred within 10 s of other probable RFI 
events, as these events tended to occur in bursts on this 
timescale. 

The effective observing time of the experiment is re¬ 
duced by the possibility of a non-RFI event being incor¬ 
rectly excluded by these cuts or by the real-time antico¬ 
incidence filter, and by the dead time of 51 ms for storing 
each detected event, during which we were unable to re¬ 
spond to further events. After allowing for these losses, 
the effective observing time was 85.5% of the total on- 
Moon time, or 127.2 hours, of which 27.8 hours lacked 
the half-limb beam shown in Fig. 1. 

After all cuts, no events remained with peak ampli¬ 
tudes in excess of 8.6<r, which is consistent with the ex¬ 
pected thermal noise background. We calibrated the flux 
sensitivity of our experiment based on the noise power in 
recorded events when pointing at the centre of the Moon, 
using a model of the lunar thermal emission [27], with a 
precision in the voltage domain of ±4.5% (random) ±4% 
(systematic). Based on this calibration, the peak ob¬ 
served amplitude is equivalent, for a pulse originating at 
the centre of a beam and detected in either polarization, 
to a threshold of 0.0047 /zV/m/MHz for the limb beams, 
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or 0.0074 /.tV/m/MHz for the half-limb beam. These 
values improve on the most sensitive previous lunar ra¬ 
dio experiment, which had a threshold of 0.0145-0.016 
HV/ m/MHz [25], 


IV. NEUTRINO APERTURE CALCULATION 

The absence of radio pulses with amplitudes above the 
thresholds given above allows us to set a limit on the 
UHE neutrino flux. We use the Monte Carlo simulation 
of James and Protheroe [28] to determine the effective 
aperture of this experiment to fluxes of UHE neutrinos, 
both diffuse and directional. In the simulation we adopt 
a lunar diameter of 32.0', its median value during our 
observations, and an Airy disk beam shape. We also 
assume constant sensitivity across the frequency range 
1.2 1.5 GHz, neglecting the minor (~ 10%) variation. 

We implement in the simulation the anticoincidence 
logic used to remove RFI, which can inappropriately re¬ 
ject lunar-origin pulses which are sufficiently intense to 
be detected directly in one beam and through a side- 
lobe of another. This anticoincidence rejection effect is 
controlled by the strictest anticoincidence cut, which ex¬ 
cludes any pulses exceeding a 4.5cr exclusion threshold in 
any beam other than the triggering beam. This effect is 
highly sensitive to the sidelobe pattern of the assumed 
Airy disk beam shape, and the receiver used in this ex¬ 
periment has been deliberately designed to minimize the 
sidelobe power [23]. We therefore scale the exclusion 
threshold by the ratio between the beam power at the 
first sidelobe of an Airy disk (1.7%) and a representative 
value for the true sidelobe power (0.5%) we derive from 
measurements of the beam pattern [29], taking the square 
root of this ratio to express it in terms of electric field 
rather than power. Note that a nanosecond-scale pulse 
detected through a sidelobe will be spread out in time, 
further reducing its peak amplitude relative to the value 
predicted by beam pattern measurements conducted with 
a continuous radio source; we neglect the minor benefit 
from this effect. 

For ease of computation, we include in the simulation 
only the radio emission from the hadronic cascade pro¬ 
duced directly by a neutrino-nucleon interaction, in ei¬ 
ther the charged-current or neutral-current case. This 
neglects the additional emission from charged leptons 
produced in charged-current interactions, which makes 
a ~ 10% contribution to the experimental neutrino aper¬ 
ture [28], smaller than the uncertainty from other effects 
[25, Apps. A & B]. Other parameters, concerning the 
neutrino-nucleon cross-section [from 30] and the radio¬ 
transparency of the lunar megaregolith, are the same as 
those taken by James et al. [25] and James et al. [31]. 


A. Small-scale lunar surface roughness 

Lunar surface roughness influences the sensitivity of 
lunar radio experiments through distinct effects at dif¬ 
ferent scales. Roughness on scales larger than the length 
of a particle cascade determines the local surface slope 
and thus how efficiently the radio emission from a cas¬ 
cade can escape the surface, which is taken into account 
by most models, including the simulations of James and 
Protheroe [28] used here. Roughness on scales smaller 
than the cascade length, but larger than the radio wave¬ 
length, causes the emission from the cascade to be scat¬ 
tered in multiple directions longitudinally along its axis, 
decreasing its intensity but increasing the solid angle into 
which it is beamed [25]. (Scattering also occurs laterally 
to the cascade axis, but this has a much reduced effect, as 
the emission beaming already has a large lateral extent.) 
The result is a decrease in the aperture for lower-energy 
neutrinos, for which the intensity of the Askaryan pulse is 
reduced below the detection threshold; but an increase in 
the aperture for higher-energy neutrinos, for which the 
Askaryan pulse remains above the detection threshold, 
and has a greater chance of being detected because of 
its increased solid angle. Both of these effects are more 
pronounced for high-frequency observations, as the wave¬ 
length is a smaller fraction of the cascade length; for an 
experiment for which the wavelength is similar to the cas¬ 
cade length, such as NuMoon [32], both effects should be 
less significant [33]. 

No complete method to determine the effects of small- 
scale surface roughness has yet been developed. A full 
model would require modelling roughness over a wide 
range of scales, including the diffractive and decoherence 
effects of propagation through the surface on the field am¬ 
plitude over a broad bandwidth, and doing so for a wide 
range of randomised surfaces. Indeed, the only existing 
estimate of the effects of small-scale surface roughness on 
lunar Askaryan pulse detection was made in the analysis 
of a previous LUNASKA experiment by James et al. [25, 
App. B], whose model we apply here. They calculate the 
aperture as a linear combination of two extreme cases, 

A(E) = rA s (E) + (1 - r)A R (E) (1) 

where the parameter r describes the relative contribu¬ 
tions from As(E), the aperture for the case of a locally- 
srnooth moon with no small-scale surface roughness, 
determined by the previous simulations of James and 
Protheroe [28]; and Ar(E), the aperture that results 
when the radio emission from the cascade is considered 
to be refracted through iVg separate surface elements 
with independent slopes. Since the surface slope will 
have some degree of correlation over the length of the 
cascade, the true aperture should lie between these two 
extremes, and r reflects the strength of that correlation. 
The effective number of surface elements N§ is taken to 
be the ratio between the cascade length and the radio 
wavelength, implying a series of surface elements along 
the length of a near-surface cascade; this parameter is 
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r = 0.56 and Ns = 12 (as above) are shown in Fig. 3, and 
corresponding limits are shown in Figs. 5 and 8. These 
show the expected behaviour, with a decrease in aper¬ 
ture at low energies and a large increase in aperture at 
high energies; the abrupt transition between these two 
regimes is an artefact of the model. However, since this 
model of the effects of small-scale surface roughness is 
less well-established than models of the neutrino aper¬ 
ture neglecting these effects, we show them in addition 
to, rather than instead of, the latter. This model will 
generally tend to overestimate the effects of small-scale 
surface roughness [25], so it is likely that the true aper¬ 
tures and limits lie between the two scenarios shown. 


V. DIFFUSE NEUTRINO SENSITIVITY 


FIG. 2. Simulation of the offset aperture ratio 
An(NsE)/As{E) in order to determine the parameter r, 
which quantifies the effects of small-scale surface roughness 
as described in the text (see Eq. 1). Points have been slightly 
offset in energy for clarity. The error bars reflect the statisti¬ 
cal uncertainty in the simulation data (inner bars) and with 
the additional uncertainty from taking Ns = 12 ± 1 (outer 
bars). The value r = 0.56 is a fit to the range E > 10 21 eV, 
for which small-scale surface roughness effects are significant. 

determined primarily by the radio frequency, with only 
a weak, logarithmic dependence on the cascade energy. 

James et al. [25, Eq. B9] show that an appropriate 
value for r can be determined by fitting the offset aper¬ 
ture ratio Ar(NsE)/A s{E): the amplitude of the elec¬ 
tric field transmitted through a single surface element 
will be reduced by a factor Ns, and the electric field is 
proportional to the primary particle energy E, hence the 
offset in energy by a factor Ns between the cases with 
and without surface elements. We use the results of their 
simulations for A-r(E) to repeat their approach, as shown 
in Fig. 2 (equivalent to their Fig. 18). We take Ns = 12, 
scaled from their value of Ns = 13 for the small dif¬ 
ference in the central radio frequency between the two 
experiments, and otherwise perform the same procedure, 
finding r = 0.56 ± 0.04. This differs slightly from the 
value of r = 0.70 ±0.06 found by James et al. [25], which 
we attribute to the significant differences between the two 
experiments (and hence between their representations in 
the simulations). 

As the simulations to determine Ar(.E') are computa¬ 
tionally expensive, we make use of the further result from 
James et al. [25, Eq. B10] that the aperture from Eq. 1 
can be approximated by 

A(E) = rA s (E) + (1 - r)k L Nl 363 As(E/N s ) (2) 

where the factor 

*l = (1 + 0.075 log^iVs)- 1 (3) 

corrects for the variation of the cascade length with the 
neutrino energy. Apertures calculated with Eq. 2 with 


The geometric apertures we calculate for a diffuse 
isotropic flux of neutrinos are shown in Fig. 3. More en¬ 
ergetic neutrinos are more likely to produce a radio pulse 
sufficiently intense to be detected in multiple beams, so 
anticoincidence rejection reduces the effective aperture 
to them. The half-limb beam improves the coverage of 
the Moon, but has a higher detection threshold than the 
limb beams, so pointing configurations that incorporate 
it show an increased aperture only at higher neutrino 
energies. 

Fig. 3 also shows neutrino apertures for other recent 
lunar radio experiments, which are based on various mod¬ 
els. The apertures for LUNASKA ATCA are, as for this 
work, based on the simulations of James and Protheroe 
[28]. The apertures for RESUN are from the analytic 
model of Gayley et al. [34], which they show to be consis¬ 
tent with these simulations. The apertures for NuMoon 
are from the simulations of Scholten et al. [35] as reported 
by Buitink et al. [32], which are more optimistic by a fac¬ 
tor ~ 10 than those determined for the same experiment 
by Jaeger et al. [36] using the aforementioned analytic 
model; James et al. [31] found a similar discrepancy rela¬ 
tive to their simulations. Nonetheless, the qualitative dis¬ 
tinctions between these experiments — the lower observ¬ 
ing frequency (113-168 MHz) of NuMoon leads to a larger 
aperture for more energetic neutrinos, while the reduced 
electric field threshold for our experiment decreases the 
minimum detectable neutrino energy — are predicted by 
all of the above models. The major remaining uncertain¬ 
ties in the apertures are associated with uncertainty in 
the neutrino-nucleon cross-section, which may alter the 
aperture by a factor of ~ 2 [25]; and the effects of small- 
scale lunar surface roughness, discussed in Sec. IV A. 

Note that the sensitivity calculations for these other 
experiments did not consider all the experimental ef¬ 
fects taken into account for this experiment, described 
in Sec. II A. Corrections for ionospheric dispersion were 
applied by LUNASKA ATCA and NuMoon, and are sub¬ 
stantially less significant for RESUN, but the effect of the 
phase of an Askaryan pulse on its amplitude has been ne¬ 
glected by all other lunar radio experiments to date, and 
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FIG. 3. Geometric neutrino aperture of this experiment be¬ 
fore inclusion of anticoincidence rejection of detected events, 
after allowing for this effect, and after also including the ef¬ 
fects of small-scale lunar surface roughness (see Sec. IV A). 
Apertures are also shown for other recent lunar radio ex¬ 
periments: a previous LUNASKA experiment with the 
ATCA [25], NuMoon [32], and RESUN [36]. See text for 
a discussion of the various models on which these apertures 
are based. 



FIG. 4. Model-independent limits to the diffuse neutrino flux 
from our experiment (solid), including the effects of anticoin¬ 
cidence rejection, and (dotted) for other experiments shown 
in Fig. 3, as well as for selected non-lunar neutrino detection 
experiments, assuming an equal ratio of neutrino flavours: 
ANITA-2 [37], RICE [38], IceCube [39], FORTE [14] and the 
Pierre Auger Observatory, combining the exposure to both 
upgoing [40] and downgoing [41] neutrinos. A range of mod¬ 
els is shown for the expected neutrino flux from GZI< inter¬ 
actions given by Allard et al. [42] (shaded), and for top-down 
UHECR origin models from Lunardini and Sabancilar [18] 
(dash-dotted) and Berezinsky et al. [17] (dot-dash-dotted). 


the signal loss due to the finite sampling rate was con¬ 
sidered only by LUNASKA ATCA. The loss of effective 
aperture due to anticoincidence rejection does not apply 
to LUNASKA ATCA and RESUN, which did not use this 
technique for RFI rejection, but does apply to the Nu¬ 
Moon experiment, which rejected coincident pulses be¬ 
tween its two on-Moon beams. Consequently, the aper¬ 
ture estimation for our experiment is more pessimistic 
than for these others, and their true sensitivity will be 
somewhat less than previously reported. 

Based on the geometric aperture A(E), we calculate 
the 90%-confidence model-independent limit on a diffuse 
isotropic neutrino flux as defined by Lehtinen et al. [14]: 

dFiso < M_ l _ (4) 

dE E t ohs A{E) y 1 

where t 0 b s is the effective observing time, and the factor 
2.3 is based on the required confidence and the expected 
Poisson distribution of the number of detected neutri¬ 
nos. The resulting limit is shown in Fig. 4 alongside 
similar limits from other neutrino detection experiments; 
and again in Fig. 5, with the effects of small-scale lunar 
surface roughness included. Note that the lunar radio 
experiments (and FORTE), which are sensitive only at 
higher energies, also observed for less time, with a maxi¬ 
mum of 200 hours for RESUN compared to 28.5 days for 
ANITA-2 and years for other non-lunar experiments. 

The spectrum of the diffuse neutrino flux depends 


strongly on the UHECR origin model. For bottom-up 
models, the UHE neutrino flux is expected to be dom¬ 
inated by neutrinos from GZK interactions. The pre¬ 
dicted flux depends strongly on the assumed source evo¬ 
lution model and on the composition of UHECRs. We 
show in Figs. 4 and 5 a range of possible fluxes for dif¬ 
ferent values of these parameters from Allard et al. [42]. 
Recent results from the Pierre Auger Observatory indi¬ 
cate that there is a trend in the composition of UHECRs 
from protons to iron nuclei at energies above 10 19 eV [43], 
implying that the neutrino flux lies towards the bottom 
of this range. 

For top-down models, there is expected to be an ad¬ 
ditional contribution to the UHE neutrino flux directly 
from the sources of UHECRs. We show in Figs. 4 and 5 a 
selection of the more optimistic models not yet excluded 
by experiment [17, 18], which predict the neutrino fluxes 
resulting from the decay of weakly-interacting scalar par¬ 
ticles called moduli, emitted from kinks and cusps in cos¬ 
mic strings. Due to the high fluxes predicted at higher 
neutrino energies, this family of models is more suited to 
being tested by lunar radio experiments. 

Since an individual kink in a cosmic string produces 
a burst of UHE neutrinos, there is a possibility that 
two or more of these neutrinos may be detected near- 
simultaneously, which would be a clear signature of this 
origin mechanism [18]. However, neutrinos which are 
too closely coincident would fail to be detected by our 
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FIG. 5. Model-independent limits to the diffuse neutrino 
flux from our experiment (solid) and others (dotted) com¬ 
pared to models for the expected neutrino flux, as in Fig. 4, 
but with the limits for lunar radio experiments incorporat¬ 
ing the effects of small-scale lunar surface roughness as de¬ 
scribed in Sec. IV A. The limit for LUNASKA ATCA is as 
reported [25], and the limit for RESUN is derived from the 
reported limit [36], which did not include small-scale surface 
roughness effects, by applying Eq. 2, adopting Ns = 13 and 
r = 0.70 as done by James et al. [25] for LUNASKA ATCA 
due to the similar detection scheme and radio frequency of 
these two experiments. For NuMoon, Ns ~ 1, and so we make 
no correction for small-scale surface roughness, displaying the 
originally-reported limit [32]. 


experiment. For two near-simultaneous neutrinos inter¬ 
acting at different points on the lunar surface, typically 
separated by ~ 1,000 km, the times of arrival of the re¬ 
sulting Askaryan pulses will typically be separated by 
> 3 ms. With this separation, the first pulse would be 
detected and identified as originating from a lunar parti¬ 
cle cascade, but the second pulse would arrive while the 
buffered data for the first event were being recorded, dur¬ 
ing which time our backend system is unable to respond 
to further triggers. If the times of arrival of the neutri¬ 
nos (and hence of the Askaryan pulses) were separated 
by > 50 ms, both pulses would be recorded, but would be 
excluded as repeated impulsive RFI. A pair of neutrinos 
could only both be detected if their separation exceeded 
the 10 s exclusion window of this anti-RFI cut. 


VI. DIRECTIONAL NEUTRINO SENSITIVITY 

The simulation we use here to determine the neutrino 
aperture also calculates its directional dependence [24]. 
The directional aperture around the Moon is shown in 
Fig. 6. We calculate the model-independent limit to 
a directional neutrino flux similarly to the diffuse case 


(Eq. 4), as 

^ dir < i_ (5) 

dE E JdtA dir {E,n,t) 

where we integrate the directional aperture A d i T (E,£l,t) 
in the direction O over the observing time. This integral, 
the directional exposure, is shown in Fig. 7. 

Centaurus A, a nearby AGN and closest FR-I type [45] 
radio galaxy, is a potential source of UHECRs and hence 
also of UHE neutrinos. As its distance [3.8 Mpc; 46] 
is less than the mean free path for GZK interactions of 
UHECR protons [47], and the neutrino flux from these 
interactions is much reduced if the UHECRs are heav¬ 
ier nuclei, the neutrino flux from Centaurus A (if it is a 
source of UHECRs) should be dominated by neutrinos 
produced by interactions in the source. Several models 
for this flux exist, determining the relation between the 
UHECR and UHE neutrino fluxes based on the assumed 
conditions in the source region. The limit that our ex¬ 
periment places on the neutrino flux from Centaurus A 
is shown in Fig. 8 along with two such models. 

Cuoco and Hannestad [50] use a semi-analytic model 
to determine the expected neutrino flux from interac¬ 
tions of UHECRs, which are assumed to have an initial 
spectrum following a broken power-law, with ambient ra¬ 
diation and matter in the core region of Centaurus A. 
Kachelriefi et al. [51] simulate the same interactions for 
a range of initial UHECR spectra, for conditions both in 
the core and in the jets of Centaurus A, although the lat¬ 
ter scenarios predict y-ray spectra that are excluded by 
observations [52]; we show their scenario for UHECRs in 
the core region with a broken power-law spectrum, which 
gives the most optimistic neutrino flux. Both models as¬ 
sume the UHECRs to be composed entirely of protons, 
and normalize the UHE neutrino flux to a UHECR flux 
based on the 2 UHECRs detected by the Pierre Auger 
Observatory originating from within 3° of Centaurus A, 
with an exposure of 9,000 km 2 sr yr [5]. In Fig. 8 we 
have rescaled the neutrino fluxes to a UHECR flux based 
on the larger-scale correlation seen in more recent data, 
with 13 UHECRs originating from within 18° of Cen¬ 
taurus A, with an exposure of 20,370 km 2 sr yr [6]. Note 
that, if these UHECRs are heavier nuclei rather than pro¬ 
tons, which the greater magnetic deflection implied by 
this larger-scale correlation supports, the neutrino flux 
will be reduced. 

The divergence at high energies between these two 
models, seen in Fig. 8, is caused by the dependence of the 
UHE neutrino flux on the maximum energy of the initial 
UHECR spectrum. Cuoco and Hannestad assume that 
there is no such cut-off, as they were concerned with the 
neutrino flux at lower energies, where this assumption 
is unimportant; while Kachelriefi et al. assume a maxi¬ 
mum UHECR energy of 10 20 eV. The true value of the 
maximum UHECR energy depends on the magnetic field 
strength in the source region of Centaurus A [53], which 
is not well constrained. 
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FIG. 6. Directional aperture around the Moon for 10 21 eV neutrinos, for a single limb beam (left) and for three beams as shown 
in Fig. 1 (right); beams are adjacent to the Moon in the directions indicated by the arrows. The effects of anticoincidence 
rejection are included, but the effects of scattering due to small-scale lunar surface roughness are not; the roughness on the scale 
of the radio wavelength has a root-mean-square value of ~ 9° [44], and so the scattering may smooth the directional aperture 
up to roughly this scale. The closest approach of Centaurus A is 28° from the Moon in the direction shown. Note that the 
maximum directional aperture is in the overlap between the regions from which a UHE neutrino could be detected by either of 
the limb beams. The exposure to Centaurus A could have been slightly increased if the limb beams were oriented to place it 
in the overlap region. 
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FIG. 7. Directional exposure of our experiment, in equatorial coordinates, for 10 21 eV neutrinos. The effects of anticoincidence 
rejection are included, but the effects of scattering due to small-scale lunar surface roughness are not. Observations were 
conducted while the Moon was on the section of its orbit marked with the thickened line, to maximize the exposure to 
Centaurus A. The ANITA experiment has substantially greater sensitivity than our experiment at this energy, but only in the 
declination range from —10° to +15°, as shown. 


Koers and Tinyakov [54] point out that if some fraction 
of the observed UHECRs are attributed to Centaurus A, 
as assumed in the above models, there are implications 
for the ratio between the directional neutrino flux from 
Centaurus A and the total diffuse neutrino flux, assum¬ 
ing that Centaurus A is a typical source of both UHECRs 
and associated neutrinos. The observed UHECR flux 
must originate from sources which, like Centaurus A, are 
within the GZK horizon, whereas the associated neutrino 


flux may originate from sources at any distance. For rea¬ 
sonable source evolution models, this implies that the 
majority of sources lie outside the GZK horizon and thus 
contribute only to the neutrino flux. Centaurus A, there¬ 
fore, should be responsible for a larger fraction of the 
UHECR flux than of the associated neutrino flux. This 
argument does not apply if Centaurus A is an exception¬ 
ally efficient source of UHE neutrinos: for example, if 
UHECRs are produced by AGN generally, but giant ra- 
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FIG. 8. Model-independent limits to the neutrino flux from 
Centaurus A from our experiment (solid), including the ef¬ 
fects of anticoincidence rejection, and (dotted) for other ex¬ 
periments from Fig. 4 that have significant sensitivity in this 
direction, including the LUNASKA ATCA lunar radio ex¬ 
periment [31]. The limits for lunar radio experiments are 
shown both with (labeled) and without the effects of small- 
scale lunar surface roughness (see Sec. IV A). The limits for 
IceCube [48] and the Pierre Auger Observatory [40, 41] are 
derived from the reported model-dependent limits for Cen¬ 
taurus A, with the energy dependence obtained for the for¬ 
mer from the aperture for the appropriate declination range, 
and for the latter from the exposure to a diffuse flux. For 
RICE [38] the model-independent diffuse flux limit has been 
scaled based on the directional dependence of the sensitivity 
provided by Besson [49], as in James et al. [31] but updated 
for the increased exposure. Models for the neutrino flux from 
Centaurus A are shown (dash-dotted) based on Cuoco and 
Hannestad [50] and Kachelriefi et al. [51]; see text for details. 

dio galaxies such as Centaurus A have a higher relative 
neutrino luminosity. 

Based on the 2 out of 27 UHECRs associated with Cen¬ 
taurus A in observational data from the Pierre Auger Ob¬ 
servatory [5], the expected UHE neutrino flux from the 
entire sky is 200-5,000 times that for Centaurus A alone, 
depending on the source evolution model, with the upper 
end of this range corresponding to the observed evolution 
of AGN [54]. Rescaling to the 13 out of 69 UHECRs as¬ 
sociated with Centaurus A on a larger scale in a more 
recent dataset [6], this range becomes 80-2,000 times. 
The exposure of this experiment to Centaurus A relative 
to the mean whole-sky exposure varies with energy, but 
is typically ~ 10 times greater. Hence, targeting Centau¬ 
rus A increased the expected number of detections for 
this experiment by l/8th in the most optimistic case, as¬ 
suming that Centaurus A is not an exceptionally efficient 
source of UHE neutrinos. 

We also consider the directional neutrino sensitivity for 
gamma-ray bursts (GRBs), which are transient sources 
of y-rays that have been proposed as sources of UHE¬ 
CRs and neutrinos [55]; the prospects for this have been 


limited by recent experimental results [56], though not 
excluded [57]. Searching for neutrinos from GRBs is a 
tempting prospect because a coincidence with a GRB 
would increase the significance of a lunar-origin radio 
pulse in our data, reducing the effective noise threshold. 
However, of the 9 GRBs in the catalogue maintained by 
the GRB Coordinates Network [58] which occurred dur¬ 
ing our observations, none were close enough to the Moon 
for our experiment to be appreciably sensitive to them. 


VII. CONCLUSION AND SUMMARY 

We have established a limit on the UHE neutrino flux 
from our lunar radio experiment with the Parkes radio 
telescope. Our experiment used a large collecting area 
and broad bandwidth, combined with a signal-processing 
strategy that came very close to the theoretical optimum 
performance (with 99.6% efficiency), giving it a lower ra¬ 
dio detection threshold than previous lunar radio exper¬ 
iments. This translates to a limit that extends to lower 
energies than previous experiments of this type, estab¬ 
lishing a tighter constraint on the diffuse neutrino flux 
for energies below 10 21 eV. Further improvement in this 
direction will require larger telescopes with broader band¬ 
width, such as the Square Kilometre Array: there is no 
substantial further improvement possible in the signal¬ 
processing performance. 

Our limit at energies above 10 22 eV, relevant for the 
detection of UHE neutrinos from some top-down models, 
does not represent an improvement over previous exper¬ 
iments. The most obvious reason for this is the incom¬ 
plete lunar coverage of the telescope beams in our ex¬ 
periment, reducing the neutrino aperture, which should 
also be remedied in future work. However, the limits for 
previous experiments have been calculated with a range 
of sensitivity models, most of which neglect several ef¬ 
fects which adversely affect the aperture. In particular, 
the possibility of a lunar-origin Askaryan pulse being de¬ 
tected in multiple beams, and hence excluded as interfer¬ 
ence by an anticoincidence filter, significantly reduces the 
aperture at high neutrino energies in our experiment, and 
probably also in others. Future analysis work for both 
past and future experiments should incorporate these ef¬ 
fects. 

There are two substantial sources of uncertainty in 
the reported neutrino flux limit. The uncertainty in 
the neutrino-nucleon cross-section in the relevant energy 
range, under standard-model assumptions, is a factor of 
~ 2, which translates almost linearly to a similar uncer¬ 
tainty in the flux limit. This cross-section is a matter of 
some theoretical interest, independent of the astrophysi- 
cal origin of UHE neutrinos, but is degenerate with the 
neutrino flux in this experiment. Our model for small- 
scale lunar surface roughness, which provides a rough 
upper bound on its effects, indicates that it improves the 
limit by a factor ~ 10 for an E~ 2 neutrino spectrum. 
Taking the uncertainty to be the range between the lim- 
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its with and without these effects, this corresponds to an 
uncertainty factor of ~ 3 around their geometric mean. 
This source of uncertainty, at least, is susceptible to fur¬ 
ther theoretical work, which should be a priority for this 
technique. 

Our observations were planned to target the nearby 
AGN Centaurus A, and achieved a sensitivity to neutri¬ 
nos from this direction which is ~ 10 x greater than the 
mean sensitivity over the sky. We have thus established 
a directional limit to the UHE neutrino flux from Cen¬ 
taurus A which is the lowest published limit for neutrino 
energies above 10 21 eV. However, this targeting makes 
only a modest contribution to the expected event rate, 
unless Centaurus A is an exceptionally efficient neutrino 
source. The targeting strategy may be better justified 
in terms of directional sensitivity to UHECRs, of which 


Centaurus A potentially contributes a larger fraction of 
the all-sky flux. 
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